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LONG-TIME BEHAVIOR OF A TWO-LAYER MODEL OF 
BAROCLINIC QUASI-GEOSTROPHIC TURBULENCE 

ASEEL FARHAT, R. LEE PANETTA, EDRISS S. TITI, AND MOHAMMED ZIANE 



Abstract. We study a viscous two-layer quasi-gcostrophic beta-plane model 
that is forced by imposition of a spatially uniform vertical shear in the east- 
ward (zonal) component of the layer flows, or equivalently a spatially uniform 
north-south temperature gradient. We prove that the model is linearly unsta- 
ble, but that non-linear solutions are bounded in time by a bound which is 
independent of the initial data and is determined only by the physical param- 
eters of the model. We further prove, using arguments first presented in the 
study of the Kuramoto-Sivashinsky equation, the existence of an absorbing ball 
in appropriate function spaces, and in fact the existence of a compact finite- 
dimensional attractor, and provide upper bounds for the fractal and Hausdorff 
dimensions of the attractor. Finally, we show the existence of an inertial mani- 
fold for the dynamical system generated by the model's solution operator. Our 
results provide rigorous justification for observations made by Panetta based 
on long-time numerical integrations of the model equations. 
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Dedicated to Professor Peter Constantin on the occasion of his 60th birthday. 

1. Introduction 

The aim of this work is to present a mathematical analysis of a highly simplified, 
or "minimal" , quasi-gcostropic model of baroclinic turbulence in a rapidly rotating 
atmosphere or ocean at large spatial scales. The model is minimal in the sense that 
it has simple but geophysically motivated representations of central factors respon- 
sible for the large-scale features of extra-tropical flows, and little else. The number 
of layers is minimal for existence of a flow-determined temperature field, and the 
model includes a spatially featureless large-scale thermal forcing, a horizontal vari- 
ation of the effective local rotation rate (through the beta-plane approximation), 
and forms of viscous dissipation at both large and small scales, without any im- 
position of flow structure by special features of the forcing or any boundary walls. 
This model is essentially the one introduced by [20] as a boundary-free modifica- 
tion of the original version of [32] , with the intent of creating a model of spatially 
"homogeneous" turbulence. However numerical integrations [30] showed the unsus- 
pected development of spatial structures (zonal jets and associated "storm tracks" ) 
of great persistence, and the presence of long time scales, and dynamically intrinsic 
spatial scales that are even now only understood by phenomenological arguments. 
Our work is intended as a step toward explaining some of the features revealed by 
the numerical integrations. 
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Although the flow in each of the two layers is two-dimensional, the dynamical 
system as a whole has certain features in common with a one-dimensional system 
that also shows spatial pattern-forming features, the Kuramoto-Sivashinksy equa- 
tion, as we explain later in this section. We then show in succeeding sections that 
arguments presented in [7], [19] and [26] for the case of the Kuramoto-Sivashinsky 
equation can be used to show the existence of a global attractor of finite fractal 
and Hausdorff dimension, as well as an inertial manifold. 

We use the version of a two-layer model that appeared in [20] and [30], which 
has two fluid layers of equal resting depth H , the lower layer has a slightly greater 
density, p2, than the upper density, pi. Motions are assumed to be quasi- geostrophic 
and to take place on a beta-plane, and the Bousinesq approximation is enforced 
(for a discussion of two- layer models see, e.g., [31]). In terms of the dimensional 
streamfunctions ^,*(x*,y*) in the layer i = 1,2, the evolution equations for the 
potential vorticity fields are: 
dQl 



g^, +Ji-^*i,Ql) = K*j,^* + i^*i-^r-^l (1.1a) 

* + Jm,Q;) = -«?.** - nljA^; + ,y*{-A)'%, (1.1b) 



dQ* 



dt 
where the potential vorticity Q* in the layer i = 1, 2 is given by: 

Qr-/o+ry* + A** + ^**, (i.ic) 

where ^* = {'^\ — \E'2)/2; and J(., .) denotes the Jacobian: 

dij) dq dtp dq 
dx dy dy dx 

K^ is the temperature (or buoyancy) damping, n'^j is the mechanical damping 
(i.e., surface drag or Ekman pumping), i/* is small scale mixing diffusion con- 
stant (or the numerical viscosity), and A is the Rossby radius defined by A = 

[gip2-pi)H/i2p2h)]'^'. ^ 

The deviations of solutions of system (1.1) from a time- invariant steady state 
with horizontally uniform vertical shear Uq has been studied numerically in [20] 
and [30]. Panetta assumed that the lower layer is at rest and the upper stream 
functions were written in the form: 

mt;x*,y*)=ri{t;x*,y*)-Uoy*, %it;x* ,y*) = r2it;x* ,y*). (1-2) 

Here ipl , ■i/'^ are the deviation stream functions and their corresponding potential 
vorticities are q* and q2, respectively. 

The transient stream-functions ipilt'jXyy) (called "eddy"), after normalization of 
ip*, q*, for i = 1,2, x* and y* using velocity and length scales Uq and A (see, e.g. 
[30]) will satisfy the following evolution equations: 

^ + J{A,qi) = -J^-iP + l)^ + ^Ti' + K-A)3V>1, (1.3a) 

^ + J{^2,q2) = -^P^\)^- ^M^i^2 - KT^ + ;^(-A)'^2, (1.3b) 

^=9(^1" "^2), (?i = AV'i-V', (?2 = Ai/'2 + V'. (1.3c) 
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Here, kt,km are the non-diniensional buoyancy and mechanical damping param- 
eters, respectively, and can have any real value. (3 is the non-dimensional central 
forcing parameter and can have any real value. The smaller /3 is in absolute value, 
the stronger the forcing. In the absence of any dissipation, non-dimensional beta 
must be less than i in absolute value, otherwise, there is no instability. Addition 
of dissipative parameters changes the stability criterion a little. 

The first two terms on the right-hand side of (1.3a) and the first term on the 
right-hand side of (1.3b) describe the interactions of the eddy fields with the im- 
posed background state; the other terms on the right-hand side in each equation 
are the dissipation terms (whether physical or artificial/numerical). We note here 
that the order of the numerical dissipation (— A)"^ is ad hoc: use of such "hypervis- 
cosity" is common in geophysical models, where it functions as a computationally 
convenient method in spectral codes to simply absorb small scale energy and en- 
strophy, and to do so primarily over a small range of wave numbers at the high 
end of the spectrum in which there is no direct physical interest. Sometimes the 
hyperviscosity is applied to the qi rather than the ^jf. from the phenomenological 
point of view, at high wavenumbers there is little difference between qi and Aipi , 
and there is no physical basis to choose one form of dissipation rather than the 
other. 

In geophysics, there are different types of multi-layer models. A derivation of 
a multi-layer model, in which the fluid consists of a finite number of homoge- 
neous layers of uniform but distinct densities, is presented in [31]. As mentioned 
above, the two-layer model is the simplest layer model that has a representation 
of baroclinic (temperature-related) dynamics important in extratropical fiows in 
planetary atmospheres and oceans. The version treated in [31] differs from system 
(1.3) in several of the linear non-conservative terms, the principal difference being 
the lack of numerical viscosity terms (i.e. v = 0). The mathematical analysis of 
multi-layer models has been studied by several authors. For example Bernier in 
[2] investigated the existence and uniqueness of solutions for a viscous multi-layer 
problem and proved the existence of absorbing sets and of a maximal attractor 
for the model; moreover, an upper bound on the dimension of the global attractor 
using the Sobolev-Lieb-Thirring lemma [34] was obtained. Bernier and Chueshov, 
later in [3] , investigated the finiteness of determining degrees of freedom for a sim- 
ilar multi- layer model. Their proof is based on the ideas introduced in [22], [23] 
and [24] (see also [14] and [15]). The models studied by Bernier and Chueshov 
arc similar to the 2D Navier-Stokcs equations and were supplemented with specific 
boundary conditions. We stress that our model is different from these models due 
to the fact that it is constructed about a background zonal shear flow, which makes 
the global nonlinear stability more involved. 

The same two- layer quasi-geostrophic model (1.3) was studied by Onica and 
Panetta in [27] and [28]. In [27], they obtained the global existence and uniqueness 
of weak solutions for tp = {ipi,ij}2) in H^^^ x H^^^ in the case of dissipative terms 
— (— A)-'^+"gi, — (— A)^+"(72: respectively, where a is an arbitrary non-negative real 
number, and k,t = 0. Later on, in [28], they continued their work and proved the 
existence and uniqueness of classical solutions that are analytic in both space and 
time. 

In the next section we will study the linear system and show that, for certain 
parameter values of the system (1.3) (mainly when the size of the domain is large 
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enough), the steady state solution -01 = ?/'2 = is hnearly unstable. However, it 
has been conjectured in [30], based on numerical evidence, that system (1.3) admits 
a finite-dimensional global attractor, in particular, that the system is non-lincarly 
globally stable, i.e. all solutions are eventually (asymptotically in time) bounded. 
This observation is an important tool in the physical and theoretical/analytical in- 
terpretations of the numerical results obtained in [30]. As we will show in section 4, 
there is no issue obtaining the global existence and uniqueness of a weak solution for 
9 = (qitQ-i) G L'^^j. X ipg^, but the main challenge is to prove Panetta's conjecture; 
i.e. proving the existence of a finite dimensional global attractor of system (1.3). 
The extra term -S^ that appears in the first equation of (1-3), which is due to the 
fact that the model is derived around a zonal shear flow, adds some difhculties in 
studying the global nonlinear stability of the system. 

The numerical results obtained in [30] suggest a stabilization effect by the nonlin- 
ear terms. Remarkably, we recall that this is a well known phenomenon for certain 
nonlinear dissipative evolution equations such as the one-dimensional Kuramoto- 
Sivashinsky equation 

dw d'^w d'^w dw 

dt dx^ dx^ dx 

and the one-dimensional Burgers-Sivashinsky equation 

dv d^v dv 

dt dx"^ dx 

subject to periodic boundary conditions, with period L > 0, in the basic period 
interval [— 7, -j], see for instance [7], [19], [26], [34] and the references therein. 

The Kuramoto-Sivashinksy equation has been studied extensively by many au- 
thors. In particular, the authors of [26] proved the global nonlinear stability, i.e. 
that the system is dissipative, in the case of odd-periodic solutions. This result 
was generalized to include all periodic solutions in [7] and [19]. See also the recent 
studies by Bronski and Gambill in [4] and by Giacomelli and Otto in [18] and by 
Otto in [29]. 

From the physical point of view, the nonlinear stability of the Kuramoto- Sivashin- 
sky equation (KSE) can be explained as follows. The linear part of the KSE is 
unstable for the low wave numbers: the low modes grow exponentially, while the 
high modes decay exponentially. The KSE without the linear part is the inviscid 
Burgers equation which preserves the energy of the solution for as long as the so- 
lution exists and is smooth. On the other hand, we know that Burgers equation 
forms a shock in finite time, which is manifested by the blow-up of the derivative 
of the solution. That is, there is a mechanism of transferring energy from the low 
wave numbers to the high wave numbers. Putting all this together, we see that the 
nonlinearity in the one-dimensional KSE pushes energy from the low wave numbers 
to the high ones, hence stabilizes the system, because the linear part is strongly 
stable for large wave numbers. 

From the mathematical point of view, the idea in the proof of the nonlinear 
stability in [26] for the case of an odd periodic solution relies on a decomposition of 
the solution w(t; x) = w(t; x) + ip, where (1 — (px) is a smooth approximation of the 
periodic delta function centered around x = ±^. The idea of translation was later 
generalized in [7] and in [19] for general periodic solutions. In space dimension 
n > 2, the problem of global existence of regular solutions and their stability 
remains open; the methods employed in the one-dimensional case do not seem to 
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be capable of extension to higher dimensional Kuramoto-Sivashinsky equations (for 
further discussion of this issue see, e.g., [1] and [5]). But, in the case of the two- 
dimensional Burger-Sivashinsky equation, the global existence of regular solutions 
follows from the maximum principle. We should mention, however, that for a 
variant of a 2D Kuramoto-Sivashinsky equation, used for the study of waves in 
fluids on an inclined plane, the author of [33] was able to generalize the methods 
developed by [7], [19] and [26], and prove that the solutions to this two-dimensional 
model are uniformly asymptotically bounded, i.e. the system has an absorbing set. 
In fact, this system has a finite-dimensional global attractor. 

Adapting similar tools and techniques developed for the one-dimensional Kuramoto- 
Sivashinsky equation, we prove the existence of an absorbing ball and a finite- 
dimensional global attractor for a slightly general version of system (1.3); that is 
when the dissipation terms arc replaced by (— A)'"V'i and (— A)™^2, respectively, 
when m > 5/2. The condition m > 5/2 is needed in our analysis in section 5 to 
prove the uniform boundcdness in time of the kinetic energy of the system. Wc 
also prove the existence of an incrtial manifold of the system and obtain an upper 
bound on the dimension of its global attractor. 

We are interested in spatial periodic perturbations about the base background 
flow with periods L in both directions. Therefore, we supplement system (1.3) with 

periodic boundary condition over the domain ^=[^7if'] , {L > 1), and initial 
data qi{x, 0) ~ qoi, i ~ 1, 2. We note that, due to the periodic boundary conditions, 
integrating (1.3a) and (1.3b) over J7 = [— ■§, -j] yield 

d f K,T f d f K^ f 

-n 11 dxdy = — / (-01 - V2) dxdy, j: I2 dxdy = — — / (^/ji - -02) dxdy. 

Furthermore, by integrating (1.3c), we also have 

gi = -- / {^i - ■1JJ2) dxdy, / 92 = :r / (^i - tp2) dxdy, 

^ Ja Ja ^ Jn 

Therefore, if we assume /„ q^ dxdy ~ 0, then J„ qi{x, t) dxdy ~ for all t > 0, for 
i = 1,2. Furthermore, the spatial average of -01 is equal to the spatial average of 
il)2 ■ Without loss of generality, we may assume that the averages over fi of both ■01 
and ■02 are zero. Moreover, we note that the system of equations: 

91 = AV^i - -(Vji - ^2), 92 = A02 + ^(^1 - 02), (1-4) 

where the unknowns are 0i and ^2 is elliptic. Specifically, if we add and subtract 
the above equations, we obtain 

91 +92 = A(0i -f -02), 91 -92 = A(0i - 02) - (01 -0'2). (1-5) 

Its clear now that system (1.5) is elliptic and thus, for given 91 and 92, we can 
solve for (0i — 02) and (0i + "02)1 under the imposed boundary condition and zero 
averages, uniquely. Thus we can solve system (1.5) uniquely for 0i and 02. The 
solutions can be computed explicitly in terms of the Fourier series. 

For simplicity of the presentation, in sections 5, 6, and 7, we will restrict ourselves 
to the case of initial data that are odd with respect to the variable y. The uniqueness 
of local solutions shows that the solution will remain odd in the y-direction for all 
time. Therefore, our solutions will satisfy the following boundary condition, which 
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is a consequence of the odd symmetry: 

X, -|) = g, [x, 1^ = 0, V'« [x, -|^ = 0, (^.T, 1^ = 0, (1.6) 

for all X e [-"I, ^] and « = 1, 2. 

The paper is organized as follows. In section 2, we study the linearized model, 
about the steady state i/'i = "02 = 0, and show that for certain range of the 
parameters this linearized model is unstable. In section 3, we introduce our notation 
and function spaces, as well as some useful regularity estimates that will be used 
later to establish our results. In section 4, we prove the global existence of solutions. 
In section 5, we prove the existence of an absorbing ball as well as the existence 
of the global attractor. In section 6, we obtain an upper bound on the fractal 
and Hausdorff dimensions of the global attractor, and we prove the existence of 
an inertial manifold in section 7. Section 8 is an appendix in which we prove a 
Lieb-Thirring type inequality. This inequality provides a useful tool in obtaining 
an upper bound for the dimension of the global attractor of system (1.3). 

2. The linear problem 

In this section, we study the stability of the linearized version of system (1.3) 
about the trivial state ipi = 1^2 = 0. We show that the linearized system is unstable 
for certain ranges of the parameters, which is an indication of a nontrivial dynamics 
of the full nonlinear system (1.3). The linearized system is given by: 

t-^-c^'t --*-'-''*■ '">' 

^ = -(/3 - ^)^ - nMAi^2 - «T0 + ^(-A)'02, (2.1b) 

^ = (01 - i'2)/2, qi = AV^i - 0, 92 = AV;2 + 0, (2.1c) 

/ q, dxdy = / ■0; dxdy = 0, i = 1, 2. (2. Id) 

Jn Jn 

We write the Fourier expansions of qj and ?/)j, j = 1,2, 



V'2k+ 2('01k-'^'2k), 



whe 


ireZ* 


2=Z2- 


(0,0), 


. From 


(2.1c), we 


obtain 








= - 


/2^|k|\ 


2 

■01k- 


-^(^1 


Lk - 02k), 


'?2k = - 


/27r|k| 


1 02k - 


qik 


[ L ) 


V L 


for: 


every 


kez*^ 


. Solving for 


07ik and 0'2k, we have for k 


1.^0, 






■01k = - 


ttk^ik 


- 7kg2k, 


02k == 


-Tk^lk 


- ak'72k 



(2.2) 



V -^ / "^/^ 1/2 , , 

"k = -J ; -2, and 7k = -j ; -^. (2.3) 



where 

'2£|k|^ ^ 
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From the above one may conclude that there exists a constant co{L, s) such that 

J2 ikp+^^(ivikp + iv2kn<co(i,s) Y. ikp^(kikp + i?2kn, (2.4) 

kez*2 kez*2 

for all s € M. It can be shown easily that system (2.1a) and (2.1b) is equivalent to 
the infinite system of ordinary differential equations: 



dq 



Ik 



dt 



= ak^ik + fefc92k, 



dq: 



2k 



dt 



= Ckgik + dkq2W, 



(2.5) 



for aU k = (fci, fc2) G Z*^, where 
27ri/ci 



Ok 



6k = 



2'Kiki 



(1 - (/? + l/2)ak) 



A€7^ 



(ak - 7k) 



Kt , . / 27r|k| 



^ (/3 + l/2)7k+^(ak-7k)-J^ 



Ck = ^^ — (P - l/2)7k + ^(ttk - 7k) - Km 
2TTiki 



L 

/27r|k 



27r|k| 
L 

7k, 

2 



ak, 



V L 



.. ^ (/3-l/2)ak-f(ak-7k)-.M(^ 

For L > 1, we note that for |k| > 1, we have 

/ L 



7k - i^ 



ai^- ly 



27r|k| 
L 

27r|k| 



7k, 



Ok- 



ak + 7k=f^3jjj , and |ak-7k|<l- 
Let Too € N, be fixed, and assume that |k| < toq, hence, we have 



(2.6) 



(ttk - 7k) 



j.|M!^] ("k + 7k) < ^ + J^(27rTOo)^ 



and 



'SA/ I — ^ I (ak + 7k) + |^(ak-7k)| + i^ I — ^ 



(ak + 7k) 
< Km + Kt + iy{2Tnno)^. 



Note that the bounds above are independent of L > 1. Hence, for any fixed 
e > 0, and for mo fixed, there exists a 6{mo) > 0, independent of L, such that if 
Km + kt + v < (5(too), then the real parts of Ok, 6k, Ck and dk can be made smaller 



than |. Therefore, if we consider the matrix 



A/k = 



flk 6k 
Ck dk 



simple calculations will show that 



tr (Mk) = flk + rfk = - 



2'Kiki 
L 



(l-2/3ak) + r?i(e). 



(2.7) 
(2.8) 



det (A/k) = Re (ok^k - Ck6k) 

= (^) [(/3-l/2)ak+(a^ + 7^)(/3'-l/4)]+r72(£), (2.9) 

with |?7i(e)| < e, i = 1, 2, ry^ are independent of L > 1, whenever km+i^t+v < (5(fco). 
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Recall that the eigenvalues of Mk are given by l/2(tr (A/k))+\/^k), where 

Ak = Re ((tr (Mk))^ - 4det (A/k)). (2.10) 

One can check that 

Ak = - (^)' ((1 - 2/3a,.)^ + 4 (ak(/3 - 1/2) + (7^ - a^)(/34 _ 1/4))) + ^3(,) 

where |r/3(e)| < e. Since ^ = ^ ('i - (^)^'), then 

'"'^^%^f(l-4/3^)-4fi-r^) 1 |+.3(e). (2.11) 



L y "'^ r '^ ' \ 2 \ L 

Now, we observe that if we choose |/3| <C 1/4, and choose i > 1, large enough, such 
that 

then, one can easily see that, under these conditions, Ak is non-negative. Therefore, 
one of the eigenvalues of Mk will have a positive real part. This proves the following 
proposition. 

Proposition 2.1. The linear system (2.1) is unstable for km,kt, and v small, 
and L large enough. 



3. Function spaces and notation 

Let n = [—^,^] ■ We denote by Lpg^(r2) the Hilbert space of fi— periodic 
functions / defined on M.'^ such that 



f\neLl^,{n) and f f dxdy = 0. 
Jn 



Let II • 11^2, (•,•) be the norm and inner product, respectively, in Lpg^(il). We 
also denote by ijpg^(r2), s e M the space of fi— periodic functions with zero mean 
characterized by their Fourier expansion: 



Hperi^) 



27r|kr ^' 



u:u=Y^ Uke^"''",Uk = M-k, Yl ( ~L ) l^'kP < 00, wq = 

We denote A = —A, with domain D{A) ~ Hp^^{Vt). When D{A) is considered as a 
subset of ipgj.(r2) with the i^ topology, the operator A : Lp^j.{il) — >• Lpg^(fi) is an 
unbounded, self-adjoint, positive operator, with compact inverse A~^. For 

«(x) = J2 ^^^'"'"^ € ^P-(^) n H;^rm, 

kez*2 
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have 



M^)= E 



ke2 



27r|k| 
L 



UkS 



2TTii 



At^x)= E ('^) ^^^"'''"^' (3-1) 



kez 



L 



for any s e R. The space ij* (fi) is a Hilbert space with the norm 



\H= 



\Aiu\\l.. 



E 



kez 



p^|k| 



Uk < oo. 



(3.2) 



We also denote EI = L^^j, x L^^j. and V'' = i/^g,. x i/p^^. We denote Ay = — ^ 
and its domain H!l^^{dy){VL) the space of functions that are characterized by: 



(2'Kk\ 



2s 



= E "fc(a;)e^" "^^ ' E ( "7~ ) l"fc(^)l^ < °^' ^01' ^11 ^ ^ 

fcgz feez ^ ^ 



2"' 2" 



For u £ H5,^j,{dy), we have 



We also define 



k£Z 



27rfc 

IT 



4u(x)=E(^) ^fc(^)e''^^' 



\\f{^,-)\\h--^ / f{x,y)dy, 



and 



ffi(ri) := < u : u e H^^j.{dy), / u{x,y) dy = 0, for all x £ 



L L 
2"' 2" 



(3.3) 
(3.4) 

(3.5) 
(3.6) 



This space is a Hilbert space with the norm ||u(a;, .)||^^ := \\Ayu{x, .)||f2. 

Remark 3.1. In this paper, C represents a positive dimensionless scale-invariant 
constant that may change from line to line. Ci represents a positive constant that 
may depend on the parameters /3, km, kt, v, m, and L. 

Let s > 0, and assume that (Qi, Q2) G V'*^^. Let (^I^i, ^2) be the unique solution 
of the elliptic system 

gi=-A*i-i(*i-*2), Q2 = -^*2 + ^(*l-*2), (3.7) 



'^idxdy ==0, i = 1, 2. 



(3.8) 



By elliptic regularity, one has (^1,^2) G V*'- Moreover, one can easily prove the 
existence of a scale-invariant positive constant C such that, for i = 1,2, and s > 
r > 0, 



|*,|k2 < CL\\A^/HA\l2 < CL^WA^ 



i\\L2, 

\A-^,\\l2<CL^^'--^\\A'^,\\l2, 



(3.9) 
(3.10) 
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and 

i(l|Avi/i||i. + i|Avi/2||i.)<IIQilli^ + IIQ2lli. 

< C(l + L^) (P*i||i. + ||A*2lli2) . (3.11) 

We recall the following Sobolcv interpolation inequality. For G H^^^{Vl), there 
exists a positive scale-invariant constant C such that 

\\A4>\\l2<C\\A''^^lA\A"^4>\\v^- (3.12) 

We also recall the following Agmon's inequality in space dimension two (see, e.g., 
[9]). For <j) G Hp^j.{rt), there exists an absolute positive constant C such that 

I!</>I1l~<c||0||^/^II0II^4'. (s.is) 

We now define the trilinear form b{-, •, •) by 

bi^,Q,Q) - J J{^,Q)Qdxdy = J (^y' ~ I^S) ^'^'"^^' ^^"^^^ 

whenever the integral makes sense. We have in particular the following property 
&(*, Q, Q) = for every * e i^per(^) and Q £ i^por(^)- (3-15) 



4. Global existence of solutions 

In this section we consider the slightly more general system of nonlinear equations 
than system (1.3) 

^ + J(^„,,) = -^-(/3+^)^ + -T^ + -^™^i (4.1a) 

^ + ./(V'2,g2) = -(/? - \f-^ - ^mM'2 - «T^ + ^^"^2, (4.1b) 

qi^-AtPi-i;, q2 = -Ai;2+i>, V' = 2 (^i " "^2), (4.1c) 

gi(-,0)=q?, q2{;0)=ql (4.1d) 

where m > 5/2, is a given positive real number, and (91,92) i^ given in H. This 
choice of m implies that the dissipation in the model is slightly stronger than 
(— A)'^"(7,;, 1 = 1,2, plus lower order terms. This system is an obvious generalization 
of system (1.3), in the sense that it contains a more general dissipative terms, which 
may or may not be symmetric. 

First, we shall say few words about the existence of solutions of the nonlinear 
system (4.1). Following ideas from the theory of the 2D Navier-Stokes equations, 
it is not difficult to show that system (4.1) is well-posed for all times (see e.g., 
[27] and [28]). However, a direct proof, which is based on energy estimates for tpi 
and Qi, i = 1,2, for arbitrary large intervals of time, does not yield the uniform 
boundedncss of solutions in time, and does not preclude a possible exponential 
growth in time of solutions. A refinement of the proof in the next section will show 
that solutions are actually uniformly bounded in time, and the dynamical system 
associated with the two-layer problem admits a compact global attractor. 
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Remark 4.1. The arguments and the estimates in this section are formal and they 
can be justified rigorously following the usual Galerkin approximation procedure 
and passing to the limit using the appropriate compactness theorem of the Aubin- 
Lions type (sec, e.g., [9], [25], [34]). 

It is clear that: 

(ktV',V'i) - (kt^,V'2) = 2kt\\iJj\\12, and {kmAtP2,^P2) = KM\\A^^^ij2\\l2. (4.2) 

Taking the L^ inner product of (4.1a) with — -01 and (4.1b) with —^2 and adding 
the equations we obtain: 

1 d 



(llAi/^^ilp + \\A'/H'2\\h + m\h) +^ E ■■^'"^'V'^lli^ + 2«T||vii2 



i=l,2 

+ km\\A'/^H\1^ = {^,^i). (4.3) 

Since by integration by parts we have (^,-01) = — (^^,gi), by the virtue of 
estimates (3.10) and (3.11) we have 
1 d 



[\\A'^'i;,\\l. + \\A'/'^2\\h+2\mh)+.Y: \\A"'"Uh 



i=l,2 

< \\A'/^ML4qi\\L- 

< C{1 + l^Y'^u^~4a^/^i\\lA\a"''^Ml^ 

< ^P-ZVilli. + g(l + ^'j^''"'V v^^i||i., (4.4) 

which can be written as 

4=1,2 

< Ci [WA'^'Hlh + U'^'Mh + ni'Wh) , (4.5) 

where Ci — Ci{v,m,L) = ^ — . Therefore, given initial data ■0° — 

ilji{0),i = 1,2, in iJpgr, and T > fixed, the solutions of system (4.1) will sat- 
isfy 

ib,eL^{0,T;H'p^,)nL\0,T;H;i,), * = 1,2. (4.6) 

Now, we observe that: 
KTii>,qi) - kt(V', 92) = KT(V',gi - 92) = -2kt('0,^'0 + V') 

= -2.CT||Ai/2^|li.-2KT||^|Ji., (4.7) 
and that 

KM{Alp2,q2) = KAf(^V'2,-^'02+V') = -KAf 11^^211^2 + K m (Atl)2 , 1p) 

<-^\\AMh + ^Mh, (4.8) 
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as well as that 

^ Y. il^^^+'^/'V'^lli^ - 2||A™/2^||i.. (4.9) 



Multiplying (4.1a) by gi and (4.1b) by g2, integrating over J7, and adding the 
equations we get: 



y J2 p("+i)/Vji2 +2j.||A™/V'!li2 +2KTpi/2^||i. 



<-(/? + :^) ( ^>9i )-(/?- ^) ( ^,-72 ) + ( ^ - 2.t) ll^lli.. (4.10) 



Adding the equation (4.3) and (4.10) implies that 



\ 2 ^ 



1 d 



otA T.M\l^+Y.\\^"'"^^\\l-^ + ^\\^\\ 



2dt 



J2 (||A("+i)/V.|li= + P^/V^.lli^) +2P'»/V'| 



<-^^-^n't-^)-w-^)m- 



^-2nT 
2 



mh 



< 



1) ^ \\A'/'A\\L4<1^\\L 



2 + 



Km 



2kt 



1^1 



2 



<^E ((l/^l + i)'ll^'/'^^lli^ + ll*lliO + 



^-2^T 

2 



ll^lli- (4.11) 



Then, 



d 



tA Ell*lli^+Ell^'^'^'''lli^ + 2IIVi 



dt 



L2 



<c,\T H\\h+J2 \\A'/''H\l2 + m'"- 



^2 I / ^ 

, i=1.2 



L2 



(4.12) 



where C2 = C2(/3,i, ka/, kt)- 

Therefore, given the initial data q^ = qi{0),i = 1,2, in L^^^., and T > fixed, 
the solutions of system (4.1) satisfy 

q,eL°°iO,T;LlJnL\0,T;H;i-'), * = 1,2. 

The estimates (4.4) and (4.12) play an essential role in proving the global existence, 
uniqueness, and the continuous dependence on initial data of solutions (see, for in- 
stance, [9], [27], [28], and [34]). In particular, this allows us to define the semigroup 
of solution operators 

S{t): (q°,q2)eHK^(gi(t),g2(t))eH. 

We summarize the above results in the following: 
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Theorem 4.2. Let (V'i,'02) ^ ^^ ^^ given. Then there exists a unique solution 
("01, "02) of system (4.1) such that 

(V'l, ^^2) e C([0, T]; ¥i) n 2.2(0, T; V"), /or aZ/ T > 0. 

Remark 4.3. We remark that the above result was established by Onica and Panetta 
in [27]. They gave a rigorous proof for the well-posedness of system (4.1) in the 
space ip = iipi,')p2) € ffper ^ Hp^j., when kt = 0- 

Theorem 4.4. Let ((j'i,(7§) G H 6e given. Then there exists a unique solution 
{qi,q2) of system (4.1) such that 

(91,92) eC([0,r];H)nL2(0,r;V™-i), forallT>0, 

and the semigroup 

S{t): {q1,q",)em^{q,{t),q2{t))em. 

is continuous from H into V™, for all t > 0. 

Remark 4.5. Further regularity results were established in [28]; the authors proved 
the Gevrey regularity of solutions following the work of Foias and Temam [17] and 
its generalization to nonlinear analytic parabolic equations in [12]. 

5. Absorbing Sets and Attractors 

In this section, we will follow the Nicolaenko, Schcurcr, and Temam (N-S-T) 
idea that was used in [26] for the one-dimensional Kuramoto-Sivashinsky equation 
in the case of odd periodic solutions, to improve the estimates introduced in the 
previous section. We notice that if the initial conditions gj and 92 ^-re periodic 
functions, odd in the variable y, then the initial conditions ipi and "02 are periodic 
functions which are also odd in the variable y. Moreover, one can easily check that 
if qi(t]x,y) and q2{t;x,y) are solutions for system (4.1), then —qi(t;x,—y) and 
—q2(t;x, —y) also satisfy the equations with the same initial values, —qi{x, —y) = 
qi{x,y) and -~q2(x,—y) = ql{x,y). By the uniqueness of the solutions of system 
(4.1) we conclude that qi(t;x;y) = —qi{t;x,—y) and q2{t]x,y) = ~q2{t]x,—y). 
Consequently, the solutions gi and 92 are odd in the variable y. A similar argument 
will show that V'l and '4'2 are also odd in the variable y. So the space of functions 
that are periodic and odd in the variable y is invariant under the solutions of system 
(4.1). Thus, we will restrict ourselves in this and the following sections to the case 
that tp = (V'1,'02) and q = (91,(72) arc periodic functions which are odd periodic in 
the variable y in fl. In this case, we will prove the existence of an absorbing ball 
(i.e., the solutions remain uniformly bounded in time) as well as the existence of a 
global attractor for system (4.1). 

Following the N-S-T idea in [26], let ipiy) be an odd C°° L-periodic function on 
n that depends only on the y variable that will be determined later. We set 

■ipi(t;x,y) = tp{y) + 'ii{t;x,y), qi{t;x,y) = q{y) + Qi{t;x,y), (5.1) 

where 

q{y) = -Ayi^iy) - ^, Qi(i; x, y) = -A*i(t; x, y) - *(t; x, y), (5.2) 

^(i; X, y) = -(*i(i; a;, y) - ^^{t; x, y)). (5.3) 
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After we substitute back, the nonlinear system (4.1) in terms of ^fi, ip2, Qi, and (72 
will become 



dQi 
dt 



+ J(*i,gi) + J(*i,g) + J(V5,Qi) 



dQi ,^ 1,9^1 ^ 



(5.4a) 

^ + J(V^2, 92) ^-^P-\)^+ ^mA4,2 - KT* + J^A"V.'2 + 92. (5.4b) 



1 



1 



V' 



gi = -KT'>P + vA'yip, g2 = --HTiJ, q = ~Ayilj- -, 



(5.4c) 



gi = -A*i-*(t;a;,y), 92 = -^^'2 + * + ^ V^ * = -(*i - V'2), (5.4d) 

Qi(0,-) = Q? = '??-9, 92(0,-) = 92°- (5.4e) 

Multiplying (5.4a) by —'^i and (5.4b) by —'02, integrating over ft, and adding the 
equations, wc obtain 



2di 



+ 2KTi|*|li2+KM||A?A2||i2 =6(V',Ql,*l)- 



\ 9a; 



(.gi,*i) - (52,01). 
(5.5) 



We set 



-^ /27rfc \ 

^ iy) = -2^ cos ( -^2/) , 

z — 1 V / 



fc=l 
with Af to be chosen later. Then we have: 

bi^p, gi, v&i) - f ^, Qi ) = J^ii' - i)^gidxdy 



A/ 



/27rfc 



2E' 

V /c=l 



y]+l\ ——Qidxdy 



For (x, y) G SI, we define 



e ^ 



E 

|fe|<A/ 

E 



9*1 

w^(a^;y) — -^(2;,y)gi(a;,y); 



9*1 
dx 



Qidxdy 



, 9*1 
9x 



Qidy 



dx. (5.6) 



(5.7) 



then 



w[x 



Wfe(x)e i~^, with Wk{x) = Y I _ e"T~yw{x,y)dy. {5. 



L 
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Notice that ^(a;,0) = Qi(a-,0) = 0, for all x G hf,f], so w{x,0) ^ for all 
X g [— ■§, t]- Thus we have for s > ^ (see section 4.1, Chapter 3 in [34]) : 



|fc|<M 



Yl ^'=(^) 

|fc|>A/ 



2s 



1/2 



\k\>M ^ 



x)\ \ J2 

\k\>M 



m 



1/2 



(5.9) 



and since Hp^^{dy){fl) is an algebra for s > 1/2, and since ^^ and Qi are odd 
periodic functions in the variable y, we have 



|A;|<A/ 



< cl2«-im(i/2-«) 



9*1 



dx 



(x,.) 



|lQi(x,.)||^.. (5.10) 



Hf, 



We may choose s ~ m — 2 > 1/2 and conclude from (5.6), using the estimates (3.10) 
and (3.11), that 



a^i 



ox 



6(^,Qi,*i)- 

< (72^2m-5jy^^(5/2-m) 



Ay = gi(x,.) 



(,,.-2) ^V|> 



^ y^ 7- 2m — 5 j\ r(5/2 — m) 



Ay ^ Ql 



L2 



A. 






!/ 



dx 



" dx 

dx 



dx 



LI 



L2 



< CL2™-4(1 + ^4)l/2^,j5/2-m||^m/2^^||2^^ (5_^^) 

Since to > 5/2, then we may choose M = M{L,m) large enough such that: 

CL^"'^'^{1 + ^4)1/2^^5/2-™ ^ ^^ (5_^2) 

and thus 

/a*i 



6(V',Qi,*i)- (^^,Qij <-||A"/2vi,j|2^. 



(5.13) 



Notice also that 



(.91, *i) + (52,^2) = -kt(V5,*) + i'(A;"Vi,*) - -kt{^,^2) 



< ^Ik'lli^ + ^ll^lli^ + '^WAfmh + 7l|A™/^*ill 



2 
L2- 



(5.14) 



As a result, (5.5) can be rewritten as 
d 



I (pi/2*i||i, + \\A'/^Mh + 2||*!ii.) + ^ (p"/2*i!ii. + P"/V2iii.) 

+ 2«;T||^||i2 < nrimh + MA^'y'^Wl.. (5.15) 
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We set 

E{t) ^ \\A^'H,\\Ut) + WA'/^MUt) + 2||#||i.(t), 
By (3.10) and (5.15) we conclude that 

d 



dt 



E{t) + C:,E{t) < 7, 



where, 



C3 = C3{iy,L,m,KT) = mm|KT, ^^2(m-i) } ' 



(5.16) 
(5.17) 

(5.18) 

7 = 7(IIV^IU^, ll^;"/'Vi||L^, ^, «t) - '^T^'lli^ + 2v\\A^I^\\l.. (5.19) 

This proves the existence of an absorbing ball for system (4.1) associated with 
periodic solutions which are odd in the y variable and initial condition ('0j,'02) G 

We will now prove the existence of an absorbing ball for system (4.1) associated 
with periodic solutions which are odd in the variable y and initial data [q^, 52) ^ '^■ 
From (4.11) and by the interpolation inequality (3.12), the estimate (3.10), and 
Young's inequality we have 



\jt I ^ 'I'^^'l^^ 



Y.\\A^f^ui^+2\mi. 



twh 



|A"/V. 



< 



i=l,2 



2P"/2^|| 



L2 



i=l,2 



2 



< ( (/32 + 1) + g(l + L^fL^'^'"' '^ ^ ^^2ra 



E ll^^^'V-: 



■»llL2- 



2KT\mi. 

^-2.^|)EpV2^.,|li. 

^ i=l,2 

(5.20) 



i=l,2 



Set 



mo= E ll'?dli.(t)+ E P'/V.lli.(t) + 2|iV'ili.(t). (5.21) 

'i=l,2 1=1,2 

Then, the inequality (5.20) together with the estimates (3.10) and (3.11) imply that 
VW + C^W{t) <C,J2 P'^V.IIi^, (5.22) 



dt 



1=1,2 



where 



C4 = 04(1', L,m) ~ min 



C(l+L4)L2(m-l)' (77,2™ J ' 



C5 = C5{v,L,m,KT,KM,l3) 



W' 



C(l+L4)2i2(m-3) 



CL 



2m 



— -2^T 



(5.23) 



(5.24) 
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Now we choose Ce = CQ{i',L,m,KT) = min{C3,C4}'- then estimate (5.17) and 
estimate (5.22) yield 

j^E{t)+CeEit)<j, (5.25) 

^W{t)+CeWit) < C5 ^ \\A'/'H\h. (5.26) 

i=l,2 

Applying Gronwall's Lemma to (5.25) implies that: 

E(t) < £;(0)e-^«* + -^(1 - e"^«*). (5.27) 

Since ipi = ip + '^i, then we have: 

\\A'/'Mmh < 2\\A'/H,ml.. + 2\\Al/'^p\\l., (5.28) 

for alH > 0. Consequently, 

WA'/'Mmh + WA'/'Mmh < 2|lAf ^"lli. + 2£;(0)e-^^* + ^(1 - e"^^*), 

(5.29) 
for aU t > 0. From (5.26) and (5.29) we have 

Applying Gronwall's Lemma yields 

W{t) < Ty(0)e~'^«* < Cs (£;(0)te-'^«* + 5(1 " e^'^"* - C6ie~<^«*) j , 
for alH > 0. Then we conclude that: 

Y. M\h{t) < M^(0)e-^«* + Cs fi?(0)te-^«* + ^{1 - e'^^* - Cgte-^^*)) , 

i=l,2 V <-6 / 

for all t > Q. Consequently, we obtain: 

hmsup J2 ll^lli^W < UC5 -: p^. (5.30) 

*-'°° j=l,2 '-'6 

This proves that the ball centered at with radius 2p in EI is an absorbing ball 
for system (4.1) when ((71,(72) ^ ^- Moreover, by the estimates (3.10) and (3.11), 
inequality (5.26) and inequality (5.29) imply 

l"--"' + c(i + J)L»...-.. Eii-^'^^-'-iii- 

< 2C5 (\\A^^^4'\\h + E{0)e'^-' + ^(1 ~ e-'='«*)') . 
After we integrate with respect to time we get: 

W{0) + 2C5 / (wA^^mil^ + 2E{0)e-^^' + ^(1 - e-^"^)") ds. 
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Thus, we conclude that 
1 



y\\\A''\^\Us)ds 




i=l,2 




^CC,(1 + L^)L^(-) /|l^,/,- 


'lli 



(5.31) 



lim sup 

t—>-oo t 

SL 

Usmg the general theory of existence of attractors (see for instance [9], and [34]) 
and the estimates above, we show the existence of a compact global attractor for 
the dynamical system associated to our model. The existence of an absorbing ball 
for the system in L'i^^(yi) together with the fact that the semigroup that generates 
the flow 

S{t): (g?,q2")eH^(gi(i),<Z2(t))eH 
is continuous from EI into H, for all t > implies the existence of a global attractor 
A. for the system in Lpg^(ri). Since our system is dissipative, i.e. the semigroup S(t^ 
is continuous from H into V"~^, and since by Rellich Lemma V™~^ is a compactly 
embedded in H, the global attractor A is compact in H. Moreover, the global 
attractor A is connected as well. (For a complete discussion of global attractors we 
refer the readers to [9], [34], and the references therein). 

Theorem 5.1. The dynamical system induced by (4.1) associated with periodic 
solutions which are odd in the variable y inft and initial data (gj, ^j) ^ H possesses 
a global attractor A which is maximal, connected, and compact in H. 

Moreover, by a similar argument as above and inequality (5.17), which implies 
the existence of an absorbing ball of system (4.1) associated with periodic solutions 
which are odd in the variable y and initial data (V'i,V'2) G V^, we conclude the 
following theorem. 

Theorem 5.2. The dynamical system induced by (4.1) associated with periodic so- 
lutions which are odd in the variable y inVL and initial data ("01: V'2) £ '^^ possesses 
a global attractor A which is maximal, connected, and compact in V^ . 

We can also prove the existence of an absorbing set in H*, s = 1,2,3,..., 
using similar methods as those in, for example ([34]), but we leave the details 
to the reader. Our next step is to obtain an upper bound on the Hausdorff and 
fractal dimensions of the global attractor A. We can also prove the existence of an 
absorbing ball in Gevrey space following [17] and the generalization presented in 
[12]. In particular, one can show that the solutions in the attractor are spatially 
analytic, and have their Fourier coefficients decay exponentially fast. 

6. Upper bound on the dimension of the attractor 

Wc now turn our attention to estimating the dimension of the global attractor 
of the dynamical system associated with (4.1). We start with the study of the 
linearized equations. System (4.1) can be written in the form: 

^+3{Rq,q)+L,q + L2q-iyA"'Rq^0, g(0) = go, (6.1) 

ot 
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where q = (<?i,'?2)^, and Rq = R{qi,q2) := -0 = {1^1,1^2)'^, is defined to be the 
unique solution of tlie coupled elliptic system 

qi^ -Aipi-ijj, q2^ -A'ip2 + 'ip, "0 = ■^(V'l - "^2), (6.2) 

ipi dxdy = 0, / V'2 dxdy = 0; (6-3) 

in Jn 

and 3{tjj, q) = (J(?/'i, qi), J(V'2, 92)), and the linear operators Li and L2 are defined 

by: 

(6.4) 
System (6.1) can be written in an abstract form as 

§=^9)9, (6.5) 

where F{q) = —3{Rq,.) — Li — L2 + vA^R. The first variation equation about 
q, where q is assumed to be the solution of system (4.1), is given by (here Q = 

(Qi,Q2)) 

^^F\q)Q, g(0)=e (6.6) 

It is easy to check that this equation can be written (formally) as 

do 

-^ + 3{Rq, Q) + J(i?Q, q) + 3{RQ, Q) + LiQ + L2Q - vA'^'RQ = 0, (6.7) 



/ dxdy = 0, Q(0) = e (6.8) 
Jn 



Following similar arguments to those in the previous sections, and the Galerkin 
approximation method, one can show the following; 

(i) if g is a solution of (4.1) with q € L°°(0,T;H) n 2.2(0, T;V"-i), for aU T > 0, 
then for ^ given in H, equation (6.7) has a unique solution Q satisfying 

g e L°°(0, T; H) n L2(0, T; V'"^i); (6.9) 

(ii) the semigroup S{t) : H ^ H that satisfies S{t)q° = q and ^^g^£. = Q{t) is 
Frechet differentiable in H, with respect to the initial values. 

We now estimate the dimension of the global attractor A, following the work of 
[8] (see also [9] and [34]). Let q = q^r) = S{T)q° be a fixed orbit (i.e. solution 
of (6.1), with time denoted by r). For fc G N, we consider ^1, . . . ,^fc, k elements 
of H, and the corresponding solutions Qi, . . . , Qk of (6.7))-(6.8), with initial data 
QjiO) = £.j, i = 1, . . . , fc. Let Pk{T) = Pfc(T, go; 67 • ■ • , 6) be the orthogonal pro- 
jector in H onto the space spanned by {(5i(t), . . . , (5a:(t)}. At a given time r, 
let Oj{t), j G N, be an orthonormal basis, such that Span{0i(T), . . . , 0fc(r)} = 
Pfc(r)H = Span{Qi(T), . . . , Qkir)}. Since Qjir) G V"-i = HL^-^ x m^~\ for a.e. 
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T, the functions 9i{t), . . . , ^^(t) also belong to V™^^, for a.e. r. We also have 

k 

Tr (F'(<z(r)) o F,(r)) = ^(F'(g(r)) o P,(r)0,-(r), 0, (r)) 
i=i 

k 

= Y,iF'{q{r))e,{r),0,{r)). (6.10) 

i=i 
Writing 0j = (0j,i,0j,2), we also define (7ju,7j,2) = ^^j, i-e., 

%i = ^7jM - 7j, ^j,2 = A7J.2 + 7j, 7j == 2 '^^J'l ~ '^^■^2^' 

/ ^j^idxdy = 0, 7j_i J7- periodic, i = 1, 2, j = 1, . . . , fc, 

Omitting temporarily the variable r, we see that 

(^'(<7)0„0,) = -{3{Rq,0j),e,) - {J{R0j,q),0j) - {3iR0„0,),0,) 

- {Li0j,0j) - {L20j,0j) + u{A"'R0,,0j). (6.11) 

Next, we note that {j{Rq,0j),0j) = {J{R0j,0j),0j) = and that 

iA^R0„0,) =(A™7,- 1, -A7,- 1 - 7,) + (A™7,- 2, -A7,, 2 + 7,) 

where 

Also, we have 

iLi9j,0j) = i-^ - '«T7j> 6'j- 1 j + i-KMAjj,2 + KT^j,0j,2) 

= (-'^t7j, -^7j,i - 7j) + (-ka/^7j,2 + KTjj,-Ajj^2 + Ij) 

= 2KT||7j||i2 + 2kt||A1/27,||2, + km|1A7,- 2||i2 - km(A7„ 2,7,), 



so, 



^(Li0„0,) = -2KT\ml^ - 2nT\\A^''-^,\\l. - km\\A^-j,2\\1^ + ^^m{A^'^1-j,2,A^'%) 



< -«:t||7,| 


h - '^T\\A'^%\\h + 


t 
\4 


"M \ 1141/2^.^1 


2 
L2- 


Furthermore, 


\iL2O„0,)\ = 


O^n)^'^-) 


+ 


(<"4)%^'^-) 





3.12) 



< (1/31 + 1)/ |Ai/2^,(r;x-,y)||0,(T;x,y)|cia;dy, (6.13) 



and 



|(J(i?0„g),0,)| 



^^^^7j (■r; a;, y)6'j- (r; a;, y) ^^/^^(t; x, y) dxdy 



< / |Ai/27,(r;a;,y)||0,(T;x,y)||Ai/2g(T;x,2/)|dxdy. (6.14) 
Jn 
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Hence, 



k 

— — fs, A -r I 1 1 Ai ' • yo ^- 1 

,4kt 



^[.Hl7,lli.+2«:T||AV27,|li.- ^--M ||AV272.,||L.J 



fe 



,=i-^n 



5] / |^i/2^,(T;x,y)||0,(r;a:,y)|(|/?| + l + |Ai/25(^;^,y)|) dxdy. 



Notice that 

^ / \A^'^^,{T-x,y)\\e,{T-x,y)\(\p\ + \ + \A^'\{T-x,y)\\ dxdy < 



3=1 -^^ 



pi/2(r; X, y)(T^'^{T- x, y) (|/3| + 1 + \A^'\{t; x, y)\^ dxdy, (6.15) 
where 

fe *; 

a(T;a;,y)=^|Ai/2^,(r;x,y)|2, p(t; x, y) = ^ |0, (r; s, y)|2. (6.16) 

i=i j=i 

By the Holder inequaUty we have 

pl/2^1/2 l^l^l ^_ I ^ |^l/2^(^. 2.^ y)|^ d^dy < 

iipiiii'ikiii^ (p'/''?iIl./3 + L'/'m + 1)) . (6.17) 

Using a version of the Licb-Thirring inequality introduced in, e.g., [9], [34], there 
is a constant C such that 

\\p\\%'<c\f:\\A'^%\\lA <C{l + Lr^'\^f2\\A^/'j,\\h] , (6.18) 
and by the version we prove in section 8, 

Ml^'<CL'^'(^f:\\A'^'l,\\h] ■ (6.19) 

Since 

\\A'/\h.,s<L'/'\\A'/\h.. 
we have 

pi/^i/^ [m + 1 + \A'/^qiT;x,y)\) dxdy < 

( " \ "" 

C(l + L4)l/4il/2^m-2 ^ ||^(m+l)/2^^.||2^ (l^/^ pl/2g|| ^, + 7,3/2^|^| ^ ^ 
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By Young's inequality we then have that 

/ pi/2ai/2 (|/3| + 1 + \A^/^q{r; x, y)|) dxdy < 

,. ™ f^(^ I r4\l/2r2m-2 , , 

^ Y: P(™+^'/^, Hi. + "^^' + ^1^ (pi/^g|li. + L^m + If) . (6.20) 

Using the version of the Lieb-Thirring inequaUty we prove in section 8, we have 



^M 

4^'^'^ 



i2\\A'/'l,M\h= f J2\^'/'l,ir;x,y)\'dxdy 

/ k \ '/' 



/f2 ^ .^1 



1/2 



,, ^ ^2r2m+2 

<-^p(™+i)/2^^.||2^ + ^2:^^ , (6.21) 



4 
i=i 

where C7 ~ Ct{kt, km)- We then conclude that 

k k ^2 T 2m+2 

Y.ip'm.0,) < - ^E ii^'"^'^''^.iii^ + ^^^ 

+ ^^ ; (Pi/'<zlli. + i'(l/?l + 1)') • (6.22) 

Now we note that, since dj, j ~ l,...,fc arc orthonornial in H, the eigenvalues 
(Aj)jgN of A™~^ satisfy 

we have (see, for example. Lemma 2.1, Chapter VI in [34]) 

k 

j=i 
for some positive absolute constant C, and therefore from the elliptic regularity 
estimate (3.11) we conclude that 



:2(m-l)"' ^ J, 



T.irm.i>,)<^j^, 



r(^ ^ 7-4U/2 r2m-2 , . 

+ ^^^^ {WA'/^Wh + L'm + ir) ■ (6.25) 
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From this we obtain 



where from the previous section, (5.31) impHes that 
C:=hmsup- / \\A^^^q{T)\\l2 dr 



(C + i2(|/3| + lf), (6.26) 



2\ 



.2^Mii^)i^,(p./V||i..±). (6.27, 

where 7 is given in (5.19), Ce = minjCa, C4}, where C3 and C4 are given in (5.18) 
and (5.23), respectively, and C5 is given in (5.24). 

Now we assume that qq belongs to the global attractor A and introduce the 
quantities rk{t) and Rk, 

1 r* 

rk{t)= sup sup {- TrF'{S{T)q„)oPk{T)dT), (6.28) 

9o6-A|Ci|<l,J=l,...fe * "'0 

i?fc=limsuprfc(t). (6.29) 

From (6.26), we have that 

,,r< r'2T2m+2 ^/i I r4\l/2r2m-2 

With the definition 

d:=TOm{fce N: i?fe <0}, (6.30) 

one can easily prove that 

"-'<( "" 'cT'" + ^" " "IT"'"" (c + I'-m . 1)-))""' < .i. 

(6.31) 

Then, from the general theory for the dimension of the attractors (see for instance 
[6], [9] and [34]), we obtain that the HausdorfF dimension of the attractor A is less 
than or equal to d. We have just proved the following: 

Theorem 6.1. Consider the dynamical system associated with the two-layer model 
(4.1) restricted to periodic Junctions which are odd in the variable y with periodic 
boundary conditions. We denote by A its global attractor, and we let d be defined to 
satisfy (6.31) where C,, 7, C3, C4 and C5 are given in (6.27), (5.19), (5.18), (5.23) 
and (5.24), respectively. C'e = min{C3,C4} and Cj = Cj^kt, km). 

^ iy) = -2 2^ cos ( — 2/ 1 > 

fc=i ^ ^ 

with M is chosen to be large enough to satisfy (5.12). Then the Hausdorff dimension 
of A is less than or equal to d and its fractal dimension is less than or equal to 2d. 
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7. Existence of an inertial manifold 

In this section wc sketch the proof of existence of an inertial manifold for the 
two-layer quasi-geostrophic model (4.1) when m > 5/2. We start by recalling the 
definition of an inertial manifold. Let iJ be a real Hilbert space and consider the 
differential equation 

ll + Mg + G((7) = 0, q{0) = qoeH. (7.1) 

Here Af is a linear self-adjoint, positive definite, densely defined operator with 
compact resolvant, while G is a nonlinear operator defined on D{M), the domain 
of the operator M. An inertial manifold for (7.1) is a Lipschitz manifold Ai that 
enjoys the following properties: (i) M is finite dimensional, (ii) M is positively 
invariant under the semi-flow generated by (7.1), (iii) A^ is a globally attracting 
set for the solutions of (7.1) at an exponential rate. For more details we refer the 
reader to the general references to the theory, and different ways of establishing 
inertial manifolds, [10] and [11]. (see also [16] and [34]). 

For the two-layer quasi-geostrophic model (4.1), we have, H = 11 = Li^^ x ip^^, 
•7= {qi,q2V e l^, and 

A/:=A"i?, G{q):=Liq + L2q + 3{Rq,q), 

where Li and L2 are linear operators and are defined by 

Li<1-[-^-'^t^,->^mA^2 + ^HJ , ^29=(^(/3 + -)^,(/3--)^ 

Recall that, due to the elliptic regularity estimate, 

\Mq\ ^ |A™-ig| (7.2) 

and thus, it is fairly easy to check that for any m > 5/2 we have 

\M-^/^G{q) - AI~^^^G{q)\ < C{K, L, m)\Mq - Mq\, (7.3) 

for all q, q satisfying \Mq\ < K, \Mq\ < K,and that G is a bounded mapping 
from D[M) into D[M^^/'^). Furthermore, if A = Xjq and A = A^v+i, we have 
\n ^ cqXiN™^^ as TV ^- 00. The spectral gap condition (see for instance [16] or 
[34]) is satisfied and the existence of an inertial manifold follows from the general 
theory (for example, see [16] and Theorem 3.1 in [34]). 

Theorem 7.1. Assume that m > 5/2, then the dynamical system associated with 
the two-layer quasi-geostrophic model (4.1) restricted to periodic functions which are 
odd in the variable y with periodic boundary conditions posses an inertial manifold. 



8. Appendix: A Lieb-Thirring type inequality 

Using a technique due to C. Foias, we will prove the following inequality (see 
also [13] and [21]). 
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Lemma 8.1. There exist an absolute positive constant C, such that for every family 
of functions 9i, ... ,9k inW^, which is orthonormal in Ell we have 

1/2 









L2 



(8.1) 



where jj = (7^.1,7^.2), J = ^t ■ ■ ik, is the unique solution of 

%, = A7,-, + .^:^(7,-i-7j-2), J jj,,dxdy = 0, i = l,2. (8.2) 

Proof. Applying the Agmon's inequality in space dimension two (3.13) we obtain 

\\A'/%\\L^<C\\A^/'ld%'\\A'^'ld]!', (8.3) 

which implies by the elliptic regularity estimate (3.11) and by the orthonormality 
oi 9i, . . . ,9k that: 

Suppose that ai,...,ak G M such that J2i=i'^'^i ^ 1- Using inequality (8.4) we 
show that 






< CL'/^ 



^a,A--7, 



l3/2. 



J = l 



1/2 



L2 



<CL'/'lj^lY^\A'/'jMy)n dxdy] 



1/4 



1/4 



thus 



Y^a,A^/^j^ix,y) 
i=i 



1/4 



(8.5) 



(8.6) 



Recall that A^^'^^j = A^^^^j^iBi + ^^^^7^,262, then the above inequality implies 
V 2 / V 2 / ^1/2 

Y^a,A^^'j,A^,y)] +(^M^/S-,2(^,y)| <CL(^p3/S|li2 
^i=i / \J=i / Vj=i 



for every ai, . . . ,ak G M that satisfy X)i-=i a? = 1- 



First, we choose 



a, = ^i/S-,i(^,y)/(E-=i(^'/SM(^,y))') 



1/2 



and later we 
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choose tti = A 



1/2 



and thus 



1/^7^,2 (.T, y)/ (e'=i (^'/'7,;2 (x, y)) ') to get that 



J2\A'^'lA^,y)\' <CL lJ2\\A'/'i,\\h 



1/2 



Y.\A'/'lA^,y)\' 



<CL|^p3/2^,| 



1/2 



(8.7) 
D 
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